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1. Introduction 

The radiative decay of J/xy — > r\ c y has been subject of extensive theoretical and experimental 
studies since many years. The current experimental results quoted by PDG [jl|] is, 

r(J/y^ r] c y) = 1.58(37) keV. (1.1) 

This value is obtained after averaging two experimental results, namely T{J /y — > T\cY) = 
1.18(33) keV by Crystal Ball [g], and the more recent value obtained by CLEOc 1.91(28)(3) keV [§]. 
The currently running KEDR experiment [Q] instead suggests a larger value, 2.2(6) keV. It is fair 
to say that the current experimental situation is unclear and dedicated charm experiment at BESIII 
is expected to clarify the situation. 

Prior to 2012 the theoretical situation concerning prediction of T(J/y — > f]cY) was not better. 
Dispersive analysis of r(T7 c — > 2y) obtained an upper bound for the width T(J/\j/ — > r\ c y) < 3.2 
keV [§]. Two different QCD sum rule calculations resulted in (1.7 + 0.4) keV g] and (2.6±0.5), 
keV An effective theoiy of non relativistic QCD found (1.5 ± 1.0) keV [§]. Lastly, two dif- 
ferent potential quark model calculatios exist, predicting an even larger value for the decay width, 
3.3 keV and 2.85 keV [|l(]]. To all of these predictions the error must be regarded not as an esti- 
mate of uncertainty intrinsic to the method, but only as susceptibility of the method to the variation 
of external parameters entering the prediction. The global picture of the theoretical predictions for 
the T(J/y — > VcY) is puzzling and inconclusive, and necessities a fully non perturbative analysis 
from the first principles of QCD. 

The first extensive study of the radiative decays of charmonia on the lattice has been reported 



in ref . [11] where the authors computed relevant matrix elements for a number of decay channels in 
the quenched approximation of QCD and with one lattice spacing only. That computation has been 
extended to the case of Nf = 2 dynamical light quark flavors at single lattice spacing in ref. [|L^]. 
In this paper we will focus on Tj c y and h c — > r\ c y, for which we compute the desired form 

factors for four lattice spacings that we could extrapolate to the continuum limit. Our result for 
T(J/y — > t~IcY) allow for a clear comparison between theory and experiment, as soon as the results 
from KEDR and BESIII become available. Our T(h c — > r\ c y) will provide us with a prediction for 
h c lifetime, and both could be compared with experimental measure when it becomes available. 

This presentation is based on our recent paper Jl3| ] where the interested reader can find details 
of all the numerically computed data entering our calculation. 



2. Hadronic Matrix Elements 

The transition matrix element responsible for the J/y — > 7] c y* decay reads, 

(r}c(k)\J™\J/ri/(p,e x )) = eJ2 c e^ vaP e* x v p a k^ 2 V(q 2 )/(mj /v + milc ) , (2.1) 

where = J2 c cy^c is the relevant piece of the electromagnetic current, with J2 C = 2/3 in units of 
e = x/4na em . Information regarding the non-perturbative QCD dynamics is encoded in the form 
factor V(q 2 ) and represents the most challenging part on the theory side. For the physical process, 
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i.e. with the photon on-shell q = 0, the decay rate is given by [11] 

8 ( A \ 

r(//V 7] c y) = — ae m {mj/y + m7]c ) ( ) \V(0)\ 2 , A = m J/v -m T}c . (2.2) 

\ m J/y/ 

Similarly h c — > rj c y transition matrix element is parametrized by two form factors Fi^(q 2 ), 
(r]c(k)\J™\hc(p,e x )) 



m \- m \ 



r 



The decay rate for the on-shell photon is Ql lQ 



T{h c -> r] c y) = 8a em (^ - ml) |F!(0)| 2 /(27mJ c ) . 



- qtl -(p + k)^2.3) 



(2.4) 



We can compute the form factor V{q 2 ) and F\(q 2 ) directly at q 1 = by using twisted boundary 
conditions [14] on one propagator, which will be labelled in the following by a superscript "0". 



3. Two-point correlation functions 

As in [15], we use the maximally twisted mass QCD [ |T6| ] gauge field configurations produced 



by ETM collaboration [17]. We extract mass of charmonia from the two point correlation functions: 



C r (t) = (£Tr [<9 r (0,0)<9 r (*,?)] ), Y {x) = c{x)Tc(x) , (3.1) 

.V 

in which the Dirac structures T are chosen as 75,^ or cTq/ , i E (1,2,3) to provide the coupling to 
the charmonium states with quantum numbers J PC = h , 1 , and l + ~, for t] c , J j\\f and h c , 
respectively. We implement the Gaussian smearing on the fermionic fields c entering |3.l| , and 



compute the quark propagators using stochastic techniques Q17|]. Charm quark mass, ;U C , at each of 
our lattices has been fixed according to the result of ref. JT~8| ] where it was shown that the charm 
quark computed from the comparison of the lattice results with the physical fully agrees with 
the value obtained by using the physical mu„ or m^. Therefore, we can say that m^ c , obtained 
by computing the effective mass from correlation function C^ c {t), and then by fitting m^(t) 
at large time separations to a constant, is merely a verification that, after a smooth continuum 
extrapolation, we indeed reproduce m e ^f' = 2.980(1) GeV. To extract the values of mj/y and m/ lc 
we proceed along the same line, computing nij t j lffh (t) using the appropriate correlation function, 
and then fitting them at sufficientely large time separation. 

In fig. [j] we show an example of two effective mass plots, as obtained by using all four lattice 
spacings explored in this work and for one value of the sea quark mass, which we choose to be 
the least light ones, for the case of v\ c and of h c mesons. We see that the effective masses for the 
pseudoscalar are excellent while the signal for h c is good but more noisy. The quality in the case 
of J/y (not shown) lies between the two illustrated. The effective masses are then combined to 

rnfit) 

m™(t) mf c (t) 
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Figure 1: Effective masses of the charmonium states, ^ (f), extracted from the two-point correlation functions at 
four lattice spacings, for one value of the sea quark mass. 
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We then fvt Rjj whc {t) plateaus to a constant Rj/ V j lc , extrapolating to the continuum limit, we get 



R J/WA - R J/y,h c t 1 + b J/yJi c m q + c J/y,hc a2 / a 3.S>] 



(3.3) 



R 



J/v 



1.0377(6) [ 1.0391(4) ] exp - 



R 



h c 



1.187(11) [ 1.1829(5) ] exp - . (3.4) 



In eq. ([O]) the parameter bj/y hc « measures the dependence on the sea quark mass, m q = 
m^ s (2 GeV), while the parameter cj/y j lc fa 3 % measures the leading discretization effects. Divi- 



sion by ag_3 9 = 0.086 fm is made for convenience. The linear fit ( 3.3 ) describes our data very well 
except for the results obtained at j8 = 3.8, that can be either excluded from the extrapolation (above 
results), or included adding a term proportional to a 4 , leading to a fully consistent result with the 
one quoted above. Having neglected the disconnected contributions to the correlation functions, the 



fact that our lattice results agree with the experimental values ( |3.4| ) can be viewed as a verification 
that they are indeed very small. 



From eq. (3.4), one can infer the hyperfine splitting: 



m J/y 



m 



I) 



(3.5) 



which after linar fit to a pametrization similar to that of eq.( 3.3) gives 

A cont - = (112 ±4) MeV [ 116.6 ± 1.2 ] PD& , (3.6) 
in good agreement with the experimental result written in brackets [flj], and in excellent agreement 



with the result of BESIII [19], A = (1 12.4 ±1.16) MeV. Note also that from the fit of our data we 
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Figure 2: Plateaus exhibited by Rj/y{t) and Rh(t) for the heaviest sea quark at /3 = 4.05. 
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find A = 1.0(3) GeV ', in qualitative agreement with ref. [20] where a tiny decrease of A is found 
while lowering the sea quark mass. Note, however, that this observation {b& ^ 0) disagrees with 
earlier findings of ref. [EFTI . 



4. Radiative Transition Form Factors 



To extract the desired hadronic matrix element (2. 1 ) we computed the three point correlators 



Qtf.t) = ^(V,H0)J] m (x)P(y)} eW-V = (£Tr [s c (y;0)tfS c (0,x)rA § (*> 



(4.1) 



where P = cy^c, V; = eye are the interpolating operators fixed at t = and t = t y = T/2 (T being 
the time extension of our lattices). Using the fact that our three-momentum q = 8/L is isotropic, 
we averaging 6 equivalent contributions, C\i, C23, C31, —C%\, — C32, — C13 to Cy, 



C v {qU) -> Jl ; V{0). 



2E Vc L 



mj/ v + m nc 



(4.2) 



where the last expression is valid for sufficiently separated operators in the correlation function ( |4. \\ . 
By renormalizing the local electromagnetic current and combining appropriately two and three 
points correlators, we can build a ratio Rj/y where we eliminate the source terms, obtaining V(0) 
from a fit to a constant as shown in fig. |2| 

Extrapolation the physical limit (wi sea = m q — > 0, a — > 0) is performed using a form similar 



to eq.(3.3). We do not observe any dependance of V(0) on the light sea quark mass. Instead the 
discretization effects are rather large, with cy = —23%. Our final result is: 

7(0) = 1.94(4). (4.3) 

We now turn to the discussion of the form factor F\ (0), relevant to the h c — > rj c y decay, as defined 
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Figure 3: Linear continuum extrapolation of the form factors V(0) and F\ (0). 
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in eq. (2.3). To that end we compute the three point correlators 

QjMt) = £(7)I(Oyr(*)PO0) &W = -<£Tr \s c (y;0) Yi YM^)ykSf(x,y)Y 5 



■v.v 



(4.4) 



Fi(0) is then obtained fitting to a constant the ratio obtained dividing the combination Cp 1 = 
[Cm +C231 +C312] /3 — [C131 +C212 + C323 +C232 + C313 + C121] /6 with the two point functions 
(see fig. ^]). We perform the continuum and chiral extrapolations, in a way analogous to eq. (3.3). 
Again, the form factor F\ (0) is insensitive to the variation of the light sea quark mass. Contrary to 
V(0), the discretization effects turn out to be smaller: we find q?, « 2%, and our final results is 



Fi(0) = -0.57(2). 



(4.5) 



5. Phenomenology 

5.1 Decays of J/y 



Concerning the radiative decay J j\\f — > r\ c y, by inserting our value (13) in eq. (2^2) we get 
r(J/\j/^>- r\cY) = 2-64(11) keV [1.58(37) keV] exf \ 



(5.1) 



where we used the measured Br(//yA ->■ T] c y) = (1.7 ±0.4)%, the full width T J/v = 92.9±2.8 keV ^ 
the physical values of m J/v = 3096.92( 1) MeV and A = 1 16.6 ± 1 .2 MeV. 

Our result for the decay rate is larger than the experimental one, and the agreement is only at 
2a. The various effective approaches presented in the introduction agree with ours too, except that 
we have smaller and controlled uncertainty. We hope more effort on the experimental side will be 
devoted to clarify the disagreement among various experiments. 
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Figure 4: Comparison of different results for T{J/y — > Tj c y). 

We note also that the quenched result of ref. [11], V(0) = 1.85(4), is only slightly lower than 
ours, while the one obtained with Nf = 2 light flavors with a single lattice spacing in ref. [12], is 
larger than ours at the same lattice spacing. 

We do not make any estimate of the size of systematic uncertainty due to the omitted s and c 
quarks in the sea. Very recently a study similar to ours has been made in ref. [ |22| ] where several 
set of gauge field configurations, obtained with HISQ action including Nf = 2 + 1 dynamical quark 
flavors. Their results are in perfect agreement with those reported here, which means that (i) the 
inclusion of the strange quark in the sea has no impact on J/y — > r\ c y and (ii) that the continuum 
results obtained by two totally different lattice regularization lead to perfectly consistent values. 

In fig. ^| we present full comparison of the experimental and theoretical findings on J/y — > r\ c y. 

5.2 h c -> 7] c y 

h c escaped the experimental detection for a long time and only recently CLEO succeeded to 



isolate this state [23] and observed that its prominent mode is precisely h c — > r\ c y, the branching 
fraction of which was later accurately measured at the BESIII experiment, with a result: Br(/j c — > 
Tj c y) = (53 ±7)% We obviously cannot compute the branching ratio on the lattice, but with 



our form factor result ( |4.5| ) we can compute the decay width using eq. ( |2.4| ). We get 

: 0.72(5) MeV. (5.2) 
-> Tf] c y) to estimate the width of h c . We obtain: 

(5.3) 



T(h c -> j] c y) 

This can be combined with the measured Br(/z, 
r(h c -> Tf] c y) 



Br(/z c -»• rjcY) 



1.37 ±0.11 ±0.18 MeV = 1.37 ±0.22 MeV, 



where the first error comes from our determination of the form factor Fi(0), and the second one 
reflects the experimental uncertainty in the branching ratio. This constitutes a prediction that will 
be interesting to check against experiment once the latter becomes available. 



To compare with other lattice results we convert the value reported in ref. [11] to our dimen- 
sionless form factor and obtain F\ (0) = —0.53(3), which agrees very well with our result. Similar 
conversion of the result of ref. [ |i2[ ] would result in Fi(0) = —0.33(1), much smaller value than 
ours, whether we compare it with the values we obtain at j8 = 4.05 or the one in the continuum. 



7 



Radiative decays of charmonia on the lattice 



Francesco Sanfilippo 



6. Summary and future perspectives 

We presented results of our analysis of the radiative decays of charmonia by means of QCD 
simulations on the lattice. Using several lattice spacings of twisted mass QCD with Nf = 2 dynam- 
ical flavors we were able to extrapolate the relevant form factors to the continuum limit. 

We emphasize that our results are obtained without inclusion of disconnected diagrams. 
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